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Abstract — For reliable and consistent quantum information 
processing carried out on a quantum network, the network 
structure must be fully known and a desired initial state 
has to be prepared on it. In this paper, we consider a 
spin network with only a single node accessible and provide 
continuous-measurement-based methods to fulfill the above 
two requirements. More specifically, the measurement results 
continuously obtained are used to update the estimate of the 
network structure, which converges to the true structure with 
high probability; for the latter problem, we further employ 
an adaptive measurement scheme to deterministically drive the 
state to a desired target state for network initialization. 

I. Introduction 

Quantum information processing can be usually performed 
on a highly networked system composed of a large number 
of subsystems [1], [2]. A specifically useful result is that 
we can achieve universal quantum computation by locally 
controlling all those subsystems [3], [4], [5], [6]. However, 
accessibility to all subsystems is impractical, in the sense 
that the cost for actuators increases with the number of sub- 
systems. Also manipulation error occurred at each subsystem 
will be accumulated, which as a result can largely degrade 
the performance of information processing. These issues 
thus stimulate development of a network control method 
based on manipulating only a small set of the subsystems, 
whereas the other subset will be controlled indirectly through 
interactions among the network. Actually some advantages 
as well as limitations of this approach in various applications 
to quantum information have been studied; for instance, 
state preparation [7], purification [8], [9], state transfer [10], 
parameter estimation [11], [12], and universal computation 
[13], [14]. Note that, when a solid system with its surface 
only accessible is employed to construct the network, the 
local accessibility is merely a constraint rather than a control 
strategy. 

Now let us turn our attention to some requirements im- 
posed on a network for quantum information processing. In 
particular, the following two are critical; (1) the dynamical 
behavior of the network must be fully known, and (2) 
a desired initial state of the network must be prepared. 
Together with the fact mentioned in the first paragraph, 
we are then reasonably motivated to develop a scheme that 
allows us to access to only a part of the network and yet to 
achieve the above two goals (1) and (2). Actually, Refs. [11], 
[12] and [8], [9] both have considered a specific network with 
limited access, and they provided a method for estimating the 
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parameters of the network and that for stabilizing a certain 
entangled state, respectively. In this paper, we especially 
consider spin networks with only a single node accessible 
and then provide continuous-measurement-based approachs 
to achieve the two goals. Continuous measurement is no 
more than the repetition of Bayesian update of the system 
state based on the measurement result, and thus it can 
suitably be applied to the problems of parameter estimation 
[15], [16], [17] and feedback control for state preparation 
[18], [19], [20]. 

We now describe our results. First, towards the goal (1) we 
particularly provide an algorithm to test whether any given 
pair of nodes of the network are connected or not; that is, 
it estimates the network structure using continuous measure- 
ment performed only on a single spin. We will demonstrate 
in numerical simulations that our scheme correctly identifies 
the network structure with high success probability. Note that 
our problem is different from that studied in [11], [12], in 
the sense that they have tackled the problem of parameter 
estimation of a certain network with known structure. We 
are rather motivated by the fact that sometimes the network 
structure itself is unknown; for instance, in the case of solid 
systems, subsystems are served by atoms produced at differ- 
ent sites possibly randomly, implying the need to estimate 
the network structure. It is also worth to note that this kind 
of structure estimation problem often appears in classical 
regime, e.g., reconstruction of the network structure of gene 
mRNA concentrations [21] and estimation of relationships 
in social networks [22], [23]. 

Next, we consider the problem of deterministic stabiliza- 
tion of the spin states polarizing along a common direction, 
i.e., \cf>, <f>, ■ ■ ■ ,<fi), where \<fi) is a desired spin state at each 
site; this spin coherent state is a most typical initial state 
required for quantum information processing in a network. 
As in the previous case, we use continuous measurement 
performed only on a single spin. As shown in [8], [9], 
any fixed measurement can stabilize a given target state 
only probabilistically; but we provide a deterministic method 
based on the so-called adaptive measurement technique [24], 
[25], which is a kind of feedback control scheme that changes 
the measured observable continuously in time, depending on 
the past measurement results. 

Notation. The dagger f denotes the Hermitian conjugate of 
a matrix. ® denotes the Kronecker product of two matrices. 
A pure quantum state is represented by a vector, which is 
conventionally written in terms of "ket" as \cf)). Especially, 
|0) = [1,0] T and |1) = [0, 1] T represent the state of spin- 
up and spin-down, respectively. Also we use |0, (f>, ■ ■ ■ ,0) = 



10)®* • • |0) for simplicity. A more general mixed state is 
represented by a non-negative Hermitian matrix p = p^ > 0; 
a pure state is a special class of a mixed state, and it can be 
expressed as p — |0)(0|. I n is the n x n identity matrix. 
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are the Pauli matrices. 

II. Network structure estimation via 

CONTINUOUS MEASUREMENT ON SINGLE NODE 

A. The estimator 

We first describe the general setup of continuous measure- 
ment. This can be physically realized by coupling an optical 
probe field to the system of interest and measuring the output 
field continuously in time. In particular when a homodyne 
detector is used for measurement, the time evolution of the 
system state p t = p\ > conditioned on the measurement 
results yt = {Y s | < s < i} is given by the following 
stochastic master equation (SME) [20], [26], [27]: 

dp t = -i[H, Pt ]dt + -fV[c]p t dt + ^H{c\p t dWt, (la) 
dY t = V7Tr[(c + c)) Pt ]dt + dW u (lb) 

where H is the system Hamiltonian, c the coupling operator 
between the system and the probe field, 7 the coupling 
strength, and dWt the standard Winner increment with mean 
zero and variance dt. Also we have defined 

T)[c]p = cpc) — - (c'cp + pc'c), 

%[c]p = cp + pS — Tr [(c + c^)p}p. 

Note that Eq. ([TJ is a quantum counterpart to the classical 
Kushner-Stratonovich equation that describes the time evolu- 
tion of a conditional probability density. Therefore, as in the 
classical case, the conditional expectation Tr (Ap) represents 
the least mean squared error estimate of an observable A. 

In our case, the system is an A-spin network whose 
structure is captured by the graph G with the set of nodes 
(vertices) V(G) and that of edges E(G); that is, each node 
represents a single spin state |0) = [a,b] T and E{G) 
represents the set of pair of spins connected with each other. 
Since each spin is of dimension 2, the whole network is 
a 2 JV -dimensional large system. Especially we assume that 
the interaction between the nodes occurs through the XY 
coupling Hamiltonian [10]: 

H= ]T A(aJ®af + aJ®^), (2) 

where crj, erj, and crj are the Pauli matrices acting on the jth 
spin; thus the notation means e.g. erj = I2 ® • • • ® o x ® • • • ® 
1%. Here for simplicity the coupling constants are assumed 
to be known and uniform. Next, suppose that only the first 
spin is accessible for manipulation; hence we are allowed to 
continuously measure only the first spin, in which case the 
coupling operator c in Eq. ([T]i is given by 

c = a\ (g> I 2 n-i. (3) 
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Fig. 1. Possible graph structures in the three spin network. In each case 
the first node is measured, as indicated by the wavy arrow. 



As a result, the state of the network with graph G, whose 
first spin is continuously measured, is subjected to the SME 
(Q]l with the Hamiltonian (0 and the coupling operator (f3). 

Now, let us consider the situation where the network 
structure G is unknown; then the problem is to esti- 
mate the network structure by determining which spins are 
connected. This is equivalent to choosing the true graph 
G"( l o) f rom a u possible nominal graphs G^,...,G^ m \ 
where m = 2 N ° 2 is the number of all combinations of 
the edges in the A-spin network. For a network com- 
posed of three spins, for instance, we have totally eight 
candidates of graphs . . . , G^ 8 ' whose edges are re- 

spectively given by E(G^) = {{1, 2}, {1, 3}, {2, 3}}, 
E(G^) = {{1, 2}, {1, 3}}, E(G^) = {{1, 2}, {2, 3}}, . . . 
and E(G^) = {}, as shown in Fig. Q] 

To solve the problem, we employ the idea developed 
in [16], [17]. Define a discrete probability distribution 
{p^, ■ ■ ■ ,2^ m '} with pW denoting the probability such that 
G^ is likely the true graph of the network. Let us further 
denote p^ the system state corresponding to the ith nominal 
graph C?W . Then, the state of the enlarged system composed 
of the classical index set and the nominal quantum systems 
is represented as 

m 

pf =J2 Pt ] \^){^\®pf\ (4) 

i=l 

where ■fJV'i)} is the set of m-dimensional orthonormal vec- 
tors with the index i corresponding to the zth graph. Now, 
the system with graph G^ is driven by the Hamiltonian 

ffW= A(6rJ®^ + cTj®^), (5) 

0-,fc)es(GW) 

while the coupling operator (f3) is commonly taken for all 
nominal graphs. That is, the enlarged system is driven by 
the Hamiltonian 

m 

H E = J2^i)(i>i\®HW (6) 

i=l 

and is measured continuously through the coupling operator 

c E = I m ®c (7) 



with c given in Eq. ©. Substituting Eqs. (|4), (O, and (0 
for the SME ([U, we have 

dp^=-i[H^,p^}dt + 7 V[c]pf ) dt 

+ ^%[c]pf ) dWi (8a) 

dp[ l) = 2^{Tr {cpf) - Tr (c E pf)}p^dWi (8b) 

dW t ' = dY t - 2^7Tr {c E pf)dt, (8c) 

where it is the measurement result; note that Y t is the output 
of the true system having the Hamiltonian H = and 
the coupling operator (0. We recursively calculate Eq. (|8} 
to update the probability distribution pW as well as the 
state pW, using the measurement result Y t ; we then expect 
that, after many iterations, p( %0 ^ will get the maximum value 
(recall that irj is the index representing the true graph). 

B. Example 1: three spin case 

Let us consider the simple network composed of three 
nodes; in this case, as depicted in Fig. [U we have m = 8 
candidates as the graph structure. For demonstration, we 
here set the true system as the chain-type network, i.e., 
qM — q(3)_ Also the initial distribution is set to the 
uniform one p^ = 1/8 Vi, since the graph structure is 
assumed to be completely unknown at the initial time t = 0. 
In this setting, we find in Fig. |2] (a) that the time evolution 
of {p t } according to Eq. (|8]l converges to the distribution 
with p( 3 ) — p^ — 0.5. That is, our scheme concludes the 
graph structure is G^ or G^ , which means that it correctly 
identifies the chain-type structure of the network. 

We now discuss about why the estimation is possible by 
measuring only a part of the network. For this purpose let 
us focus on the estimate (conditional expectation) of the z- 
component of the measured spin. In general, the continuous 
measurement of <j z tends to increase the absolute value of 
z-component of the spin state [19], while the z-component 
of the first spin is diffused among the network due to 
the XY coupling Hamiltonian [28], [29]; i.e. spin diffusion 
occurs. Hence, intuitively, if the network is relatively "small" 
and the spin wave gets back to the measured spin quickly, 
the estimate of the z-component of the measured spin will 
change very fast, while in the opposite case the estimate 
will change slowly. We plot the time evolution of = 
Tr (cpf o) ) in Fig. [2] (b) and those of zf = Tr {cpf) in 
Fig. |2] (c). These figures may support the validity of the 
above intuitive observation; that is, the chain structure is 
relatively a "large" network, and the corresponding estimates 
Z^ and Z^ are changing slowly. Moreover, Figs.|2](b) and 
(c) indicate that the time evolutions of the nominal estimates 
zj. 3 ^ and zf^ are similar to that of the true one z[ l °\ 
while the others are not. These facts clearly mean that the 
measurement even only on a part of the network certainly 
brings useful information for estimating the structure. 

Remark: The algorithm © identifies the true network 
structure with high probability (in this case, it is over 80%), 
but it sometimes fails. For reliable estimation, therefore, 
we need to repeat the procedure and evaluate the averaged 
probability distribution {(p[^)}. 
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Fig. 2. Numerical simulation of Eq. (§} for the three-spin network with 
parameters A = 7. The true graph structure is set to G^ ! °' = . Figure 
(a) shows the time evolution of {p+ }, which converges to p(^) = p( 5 ) = 
0.5. Figures (b) and (c) show the conditional expectation of c of the true 
system with graph G^°^ = and those of all the nominal systems 

driven by Eq. (§}, respectively. 



C. Example 2: four spin case 

We next consider the network composed of four spins, in 
which case there are totally to = 64 candidates for the graph 
structure. First, let us assume that the true graph is a chain 
depicted in Fig. [3] (a); the index of this graph is iq = 27. 
Figure [3] (a) shows the time evolution of (p^) averaged over 
ten sample paths. The initial distribution is set to the uniform 
one p^ = 1/64 Vi, because of the same reason explained 
in Sec. II-B. It is seen from the figure that (pf 1 ) takes the 
maximum value at i = 27,29,41,45,50,51, all of which 
correspond to the chain form. Hence, the estimator certainly 
identifies the true graph structure. 

Figures [3] (b) and (c) illustrate the cases where the true 
graph structure is given by the star and the square, respec- 
tively; in both cases, the true structure is correctly identified. 
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Fig. 3. Time evolution of the probability distribution (p| ) averaged over 
ten sample paths. The true graph structure is given by (a) chain (b) star, 
and (c) square. The parameter is set to A = 7. 



III. Network initialization via continuous 

MEASUREMENT ON SINGLE NODE 



Once the network structure has been identified, we can 
move forward to the stage of an appropriate state preparation. 
In this section, under the assumption that the network dynam- 
ics including its structure is completely known, we provide a 
scheme that deterministically stabilizes a spin coherent state 
\<j>,<f>, ■ ■ ■ ,4>) on the network, via continuous measurement 
performed only on a single spin. The deterministic driving of 
the state is brought from the mechanism of adaptive measure- 
ment; to make the idea clear, in this section, we first present 
an adaptive measurement method for state preparation of a 
single spin. Then it is applied to quantum networks with only 
a single node accessible. 



A. Adaptive measurement for single spin state preparation 

As indicated from Eq. (|T), measuring a quantum system 
always brings a stochastic driving of the state. This means 
that only a fixed measurement does not help to stabilize a 
state deterministically at a certain target state. Combining 
measurement with feedback control is thus expected to over- 
come this issue, as actually demonstrated in several studies 
[18], [19], [20]. Adaptive measurement is a kind of feedback 
control, which does not introduce an additional actuator 
for control but rather changes the detector configuration 
based on the past measurement results. Thus a merit of 
adaptive measurement may appear in a practical situation 
where it costs cheaper than adding an actuator for feedback 
control. The applicability of this scheme for state preparation 
has been demonstrated in for instance [24], [25]; we here 
consider the same problem studied in these references and 
add a new result. 

The problem of stabilizing a spin state via adaptive 
measurement is described as follows. Recall that the state 
under continuous measurement evolves in time according to 
Eq. (HJ. We then set H = and assume that the coupling 
operator c can be changed in time as a function of p t , so that 
p t will deterministically converge to a desired target state. 
Without loss of generality, the target state can be set to |0). 
Then, it is reasonable to parameterize the coupling operator 
in the following form: 

f cos 9 t 
e lSt sin0 t 



ct = 



1 sin 
i lSt sin t — cos 9 t 



(9) 



That is, the problem is to determine the dynamics of the 
parameters (0 t , St) so that the state p t governed by the SME 
(Q} with H = and Eq. © will converge to |0). 

The solution is obtained by introducing an appropriate cost 
function, as often done in control theory. Here we define the 
cost as 

J t = 1 - Tr (a z Pt ). 

This is non-negative and takes the minimum value only 
when p t — |0)(0|. Therefore, the parameters (9 t ,S t ) should 
be chosen so that the cost will decrease in the mean sense; 
that is, if d£[J t ]/dt < is satisfied for any p t 7^ |0)(0|, the 
cost decreases until it reaches the minimum value 0, implying 
that p t — > |0)(0| almost surely. For more detailed discussion 
about this stochastic convergence, see [30]. To make a control 
law, let us parameterize the state as 

p-iPt 

Pt 



1 
2 



1 + r t cos a t 
r t e t ^ t sina t 



r t e Lpt sin a t 
1 — rt cos a t 



(10) 



Then, the derivative of E[J t ] is given by 
dE[J t ] 



dt 



2 



(cos(20 t ) — 1) cosa t 
+ sin(20 t )sin a t cos(<5 t -fa) . (11) 
Hence, when choosing the tuning parameters as 

(6 t ,6 t ) = (a t /2,p t ), (-atAft+Tr), (12) 



we have 



dE[J t ] 
dt 



n 



(1 — cos off) < 0. 



Ul 0.5 



Then, from the result of stochastic stability [30], 
dE[Jt]/dt — > 0, thus equivalently a t —> or r t — > 0, is 
guaranteed. This means that after long time limit the state 
will be found on the positive half of the z axis in the 
Bloch sphere. But it is well known that an ideal continuous 
measurement of an observable always increases the purity 
of the conditional state; i.e., r t — > 1. Combining these two 
results, we can conclude that the state converges to the 
target state |0) almost surely. The adaptive law ([T2| > has been 
found in [25], though without rigorous proof. Hence here we 
present the result as a new contribution. 

Theorem 1: The spin state subjected to the SME ([B with 
H = and the adaptive measurement law (0, dTOb . and ([T2l 
converges to the target state |0) almost surely. 

B. Network initialization via adaptive measurement on sin- 
gle node 

Now we apply the adaptive measurement scheme devel- 
oped in the previous subsection to TV-spin quantum networks. 
We again consider the case where we are allowed to contin- 
uously measure only a single spin of the network. Also the 
network is driven by the XY coupling Hamiltonian as before. 
In this setting, the state under continuous measurement obeys 
the SME (Q~|i with the following system matrices: 




tb) 



ti, 0.5 
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t, 0.5 
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Fig. 4. Ten sample paths of the fidelity F t = (0 W4 \pt lO^ 4 } for several 
network structures: (a) chain, (b,d) star, and (c,e) square. In (d,e) a local 
control indicated by the red arrow, which breaks permutation symmetry of 
the network, is introduced. The parameters are set to A = 7. 



The point here is that Eq. ( fTol l has the same form as that for 
the single spin case, Eq. (fTTb . Therefore, the adaptive law 



%,8' t ) = {a! t /2,p' t ), (-04/2,$ + *) 



(18) 



H : 



(j,k)eE(G) 

( cos 6' t e -w *siii0{ 
V e<sm.6' t -cos0' t 



5V-i. (14) 



The target state can be set to the spin coherent state |0 8Ar ) = 
1 , 0, ■ ■ • , 0) without loss of generality, as in the single spin 
case. That is, the problem is to stabilize lO 18 ^) in the network 
by applying continuous adaptive measurement performed 
only on the first spin. 

Let us now determine the adaptive law of the parameters 
(0' t ,6' t ) of the measured observable (TBI . For this purpose, 
similar to the single spin case, we consider the following 
cost function: 



leads to dE.[J' t ]/dt — > as before, which thus concludes 
p' t — > |0)(0|. Also note that the measured observable (fT4T > 
+ ctJ <S> Ofc), (13) then becomes Eq. ©. 

Unfortunately, the above result does not necessarily mean 
the convergence of the whole state p t to the target. Nonethe- 
less, now we know that the z-component of the first spin, 
which is continuously raised via the adaptive measurement, 
is diffused among the network due to the XY coupling 
Hamiltonian; hence the z-components of all spins will also 
increase. Furthermore, the target state |0® ) is a steady state 
of the SME with ( [T3l and (0 (see Appendix). From these 
two facts, we expect that the state will converge to the target, 
if it is the unique steady state of the controlled SME, i.e., 
the SME driven by the adaptive measurement. 



Ji = N 



Tr (£f =1 <r» 



(15) 



which is non-negative and takes the minimum value only 
when p f = \0® N )(0® N \. The time derivative of E[J' t ] is 
given by 



dE[J' t ] 
dt 



(cos(20j) - l)cosc/ f 

+ sin(26»0 sine*; cos($ - #)] , (16) 

where (r' t , a' t , f3' t ) are the parameters of the reduced quantum 
state p' t — Tr (2,3,4,. ..,N)[Pt] a s follows: 



Pt = 



1 ( 1 



r; cos a, 



sin a f 



sin a t 



r; cos a. 



(17) 



C. Numerical simulations and permutation symmetry 

First, let us consider a four qubit network with chain graph 
whose first spin is continuously measured with the adaptive 
law ( [T8] ). A notable feature of this network is that in this 
case |0® 4 ) is the unique steady state of the controlled SME, 
hence as mentioned above the deterministic stabilization 
would be expected. To evaluate the performance, we use 
the fidelity F t = (b m \p t \0 m ), which takes the maximum 
value 1 only when p t = |0® 4 )(0 (gl4 |. Fig.ffl(a) displays ten 
sample paths of Ft, with the initial state po = (I2/2)® 4 , and 
show that all trajectories converge to 1. Actually the same 
behavior was always observed in 300 simulations. Therefore 
it seems conclusive that the deterministic state stabilization 
is attained. 

In the above example, a crucial property to realize the de- 
terministic convergence is that the target is the unique steady 



state of the controlled SME. Hence, we should study systems 
that do not have such uniqueness property. In particular, 
systems having permutation symmetry [14] are important; 
this property means that the Hamiltonian is invariant under 
the exchange of some specific pairs of spins. For example, 
the "star" formed network with Hamiltonian 

H = \(XXII+IXXI+IXIX+YYII+IYYI+IYIY), 

(19) 

where e.g. XXII = a x <g> a x <g> I2 <8> I2, is permutation 
symmetry with respect to the exchange of the third and fourth 
spins. In this case, in addition to |0® 4 ), the "entangled state" 

|Ent} = |0,0>®i(|0,l>-|l,0)) 

satisfies ff|Ent) = and is a steady state of the SME due to 
the reason given in Appendix. |Ent) is a superposition of two 
states |0, 0, 0, 1) and |0, 0, 1, 0), which is invariant under the 
permutation operation. Now the state could go toward |Ent), 
hence the deterministic convergence to the target would not 
be expected. Actually, as seen in Fig. [4] (b), for the network 
with the star graph Hamiltonian (fl9t , the state does not 
converge to the target deterministically; it converges to a 
mixed state on the subspace spanned by the steady states of 
the controlled SME. A similar result is observed in Fig. 0(c) 
where the network with square-graph Hamiltonian is studied. 

A reasonable method that circumvents the above issue 
would be the one adding a local control Hamiltonian such 
as Hioc = IIIZ, which breaks the permutation symmetry 
structure. In fact, with this local control Hamiltonian, we 
can attain the deterministic stabilization of the state at the 
target, both in the networks with star and square graphs, as 
shown in Fig. [4] (d) and (e), respectively. 

IV. Conclusion 

In this paper, we have provided particular continuous- 
measurement-based schematics for estimating and initializ- 
ing spin networks that allow access only to a single node. As 
mentioned in Introduction, these are essential requirements 
for performing quantum information processing; in fact, if 
the structure of the network is identified and an appropriate 
initial state is prepared on it, as shown in [13], universal 
quantum computation on that network is possible via only a 
local control on a single spin. 

APPENDIX 

It is directly proven from [31], [32] that a pure state |</>) 
is a steady state of the SME ( [Tal l if and only if \<f>) is a 
common eigenvector of iH + c^c/2 and c. Clearly, |0® ) 
is an eigenvector of Eq. (0. Also noting the relation (a x ® 
ct x +<jy ® erf )|00) = 0, we can readily see H\Q® N ) = 0. 
Therefore, |0® w ) is a common eigenvector of iH + 
and c, hence it is a steady state of the SME. Note that \1® N ) 
is also a steady state, but the state driven by the adaptive 
measurement always converges to a state with its first node 
given by |0). 
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